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, We present a new formula for the rate at which cosmic strings lose energy into gravitational radiation, 

5^ , valid for all piecewise-linear loops of infinitely-thin cosmic string . At any time, such a loop is composed 

bp. 

of N straight segments, each of which has constant velocity. Any cosmic string loop can be arbitrarily-well 
approximated by a piecewise-linear loop with TV sufficiently large. The formula is a sum of 0{N^) polynomial 
and log terms, and is exact when the effects of gravitational back-reaction are neglected. For a given loop, 
the large number of terms makes evaluation "by hand" impractical, but a computer or symbolic manipulator 
yields accurate results. The formula is more accurate and convenient than previous methods for finding the 
gravitational radiation rate, which require numerical evaluation of a four-dimensional integral for each term 
in an infinite sum. It also avoids the need to estimate the contribution from the tail of the infinite sum. The 
formula has been tested against all previously published radiation rates for different loop configurations. In 
the cases where discrepancies were found, they were due to numerical errors in the published work. We have 
isolated and corrected the errors in these cases. To assist future work in this area, a small catalog of results 
for some simple loop shapes is provided. 
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INTRODUCTION 

Cosmic strings are one-dimensional topological defects which appear in some gauge theories of the 
fundamental interactions. Strings would appear at phase transitions where symmetries of the fundamental 
interactions are spontaneously broken [1,2,3]. It is thought that cosmic strings might have formed as the 
universe expanded and cooled during the past. They are remarkably simple objects, characterized by a single 
parameter fi, which is their mass-per-unit length. For strings of cosmological interest, the expected value 
of the dimensionless parameter G^/c? is of order 10~^, where G is Newton's gravitational constant and c 
is the speed of light. The strings of interest for this work are strings without ends - thus they are always 
topologically in the form of circles, or possibly infinite in length (in a spatially infinite universe). 

The dynamics of a network of cosmic strings in an expanding universe have been thoroughly studied 
[4,5,6]. To describe these dynamics, it is useful to divide the strings, for the purposes of labeling, into two 
categories: the long string (length greater than the horizon length) and the loops (all the rest). Early work 
on cosmic strings established that the energy-density of the long strings was a small constant fraction (of 
order G/x/c^) of the energy-density of the cosmological fluid. In the literature, this is referred to as "scaling" 
behavior. The long string network maintains scaling behavior by constantly "chopping off' loops of cosmic 
string. This process takes place whenever long strings meet each other, make contact, and "intercommute" . 
Typically, after a loop is chopped off it begins to oscillate due to its own tension, undergoing a process of 
self-intersection (fragmentation) and eventually creating a family of non-self-intersecting oscillating loops. In 
the absence of gravitational radiation, these loops would survive forever, oscillating periodically, and would 
eventually come to dominate the energy-density of the universe [1] . However, these loops gradually decay 
away due to the emission of gravitational radiation [3] . 

The emission of gravitational radiation is thus of fundamental importance to the topic of cosmic strings. 
Indeed, the resulting stochastic background of gravitational radiation left behind from the families of small 
string loops provides the main cosmological constraints on cosmic strings, through two observable effects [7 
and references therein]. First, because gravitational radiation contributes to the energy-density, it affects 
the expansion rate of the universe. The amount of gravitational radiation must not be too great or it would 
interfere with the highly-successful standard model of nucleosynthesis. Secondly, the amount of gravitational 
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radiation must not be too great to interfere with the extremely small timing residuals observed in the periods 
of a number of carefully observed fast pulsars. The work on these cosmological constraints is reviewed and 
updated in [7]. 

During the past fifteen years a number of detailed calculations have been carried out to determine the 
rate at which cosmic string loops convert their energy into gravitational radiation. The power radiated by a 
given loop is 

^=^=7Gm'c, (1.1) 

where E is the energy radiated in gravitational waves in a single oscillation of the loop, At is the period 
of that oscillation and 7 is a dimensionless constant that depends only on the shape of the loop and its 
velocity at any fixed instant in time. Thus the problem is to determine the numerical value of 7 for a 
given string loop. Because loops are relativistic objects which have typical velocities of order c, the simplest 
approximation formulae like the quadrupole approximation are not of much use, although in some cases 
they are reasonably accurate [8]. Vachaspati and Vilenkin [9] carried out the first detailed calculation 
of 7 for a simple generalization of the circular loop. While some of the integrations were carried out 
analytically, the final integration over directions could only be done numerically. The next work was a hybrid 
analytic/numerical calculation by Burden [10], for a set of loops which were a variation of the Vachaspati 
and Vilenkin family. The first entirely exact analytic calculations were done by Garfinkle and Vachaspati 
[11], who considered a special family of "kinky" string trajectories. These are the simplest piecewise linear 
loops for which the exact formulae given in this paper may be applied directly. Additional work by Durrer 
[8] repeated some of the earlier calculations of the previous three groups and also investigated the accuracy 
of the quadrupole approximation for determining 7. The next work was a pair of papers by Scherrer, 
Quashnock, Spergel and Press [12], and by Quashnock and Spergel [13], which developed numerical and 
analytic techniques to study the effects of gravitational back-reaction on the shape and motion of the cosmic 
string loops. This is the first work which examines the way in which the shape of a string loop is changed as 
a result of the emission of gravitational radiation. (In our paper, these effects are not taken into account - we 
assume periodic motion of the loop). In addition to verifying some of Burden's results, they also obtained 
interesting results concerning the distribution of 7 for typical families of non-self-intersecting string loops. 

3 



Recent work by Allen and Shellard [14] used FFT (Fast Fourier Transform) methods to determine values of 
7 for the loops produced in their numerical simulation of cosmic string networks in an expanding universe. 

These investigations are important for the reason mentioned previously; the cosmological consequences 
of cosmic strings are largely visible via the direct and indirect effects of the gravitational waves produced 
by the string loops. Thus, "typical" or expected values of 7 appear in expressions for observable quantities 
such as the present-day energy-density expected in gravitational waves. Much of the research work on 
gravitational radiation by cosmic string loops has been motivated by a desire to determine the "typical" or 
"expected" values of 7. Thus, Scherrer, Quashnock, Spcrgel and Press [12] give a histogram of the expected 
values of 7; the mean is 7 = 61.7 and the median is 55.4. 

In much of the literature on this topic, the method used to determine 7 is numerical. A loop of cosmic 
string radiates at discrete frequencies corresponding to the different normal modes of motion, so 7 = X^^i 7n 
is a sum of terms arising from each of these normal modes, labeled by n = 1, 2, 3, • • •. The value of each 7„ 
is given by an integral over the two-sphere of a particular function. This function, in turn, is a product of 
integral transforms over the world-sheet of the loop. Except in certain highly- symmetric cases, numerical 
methods must be used to determine the required four-dimensional integrals. Because it is only practical to 
determine 7„ up to n of a few hundred or thousand, one must extrapolate the dependence on n in order 
to estimate the "tail" terms arising in the infinite sum over n. This process is error-prone because the sum 
over n may converge very slowly (if at all - with back-reaction neglected, 7 may be infinite!). Also, if the 
integration over the two-sphere is not done accurately enough, the 7„ will be inaccurate for large n. This 
will cause the sum over n to converge at the wrong rate. Indeed, we have found that many of the previously- 
determined values of 7 given in the literature are incorrect (typically by a factor of order two) because this 
tail has either not been included, or has been incorrectly estimated. 

In this paper, we develop a new method for determining 7. Our method yields an exact analytic formula 
for 7, valid for any piecewise linear cosmic string loop with piecewise linear velocity. (Equivalently, both the 
left- and right-moving trajectories are piecewise linear). This piecewise linear requirement is really not very 
restrictive, since in practice any cosmic string loop can be arbitrarily closely approximated by a piecewise 
linear cosmic string. Thus one can use this formula to determine 7 to arbitrary precision for any cosmic 
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string loop. Remarkably, our formula involves nothing more complicated than log and arctangent functions. 
However it is the sum of order N'^ terms, where N is the number of piecewise linear segments, and thus in 
practice is extremely cumbersome to evaluate without the assistance of a computer or symbolic manipulator. 
We stress that although our formula will probably never be evaluated without the use of a computer, it is not 
a numerical method, but rather is an exact formula. It is also fairly rapid - with A'' of 32 a DEC ALPHA class 
workstation can evaluate 7 in less than 2.5 seconds. We are making our computer code, which provides one 
implementation of this algorithm, publicly available via anonymous FTP from the directory pub/pcasper at 
alphal .csd.uwm.edu. 

In order to test our new formula, we constructed piecewise linear approximations to the smooth cosmic 
string loops studied in earlier published calculations of 7. In a number of cases we obtained very close 
agreement between the value of 7 given by our formula and the published values. However there were also 
a number of cases in which the results did not agree. Section 7 contains further details of these cases. In 
every case where we had found disagreement, we were able to show that our formula in fact had given the 
correct result. The disagreement in each case was due to numerical errors in the original work. Many of 
the published values of 7 are off by about a factor of two. For example, Vachaspati and Vilenkin give the 
value 7 = 54.0 for the case a = 0.5 and (p = O.Btt in equation (2.24) of reference [9]. The correct value is 
7 = 97.2 ± 2. Note that the value of 7 = 97.2 is exact (to three significant figures) for the piecewise loop 
which we used to approximate the smooth Vachaspati and Vilenkin loop. The error bar of ±2 in 7 arises 
because our piecewise approximation had only N = 64 segments. 

We intend to use this exact formula in future work, for example to identify the shape of a cosmic string 
loop with the smallest value of 7, and to repeat some of the work of Scherrer, Quashnock, Spergel and Press 
[12] concerning the distribution of 7 values of non-self-intersecting loops. 

The remainder of the paper is organized as follows. Section 2 describes the periodic motion of a cosmic 
string loop oscillating in flat space-time. It establishes notational conventions and a number of basic results. 
Note that in our approximation, the back-reaction of gravity on the string loop is neglected, so that space- 
time remains flat. In this context a given string loop oscillates periodically and radiates forever. Section 
3 starts with a standard result [15] for the energy radiated by gravitational waves emitted from a periodic 
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source, and obtains an integral representation for 7 in terms of the gravitational interaction of the cosmic 
string world-sheet with itself. This was motivated by (and is almost identical to) a calculation given in 
appendix B of [13]. In Section 4 we restrict our attention to the special case of piecewise linear loops, 
and establish notational conventions for such loops. The "corners" of the piecewise linear loop trajectory 
may be discretely labeled; their positions/velocities contain all information about the loop. The integral 
representation for 7 is then expressed in terms of these discrete quantities. In Section 5, the formula for 7 
is simplified and expressed as a sum of elementary integrals. These integrals are three-dimensional volume 
integrals; the integrand is a Dirac delta function of a quadratic form in x,y and z. These integrals are 
evaluated in closed form in Section 6. This section contains the main result of the paper, which is an exact 
closed-form expression for 7 in the piecewise linear case. Section 7 contains the results of our investigation 
of the existing literature, reporting both on those cases where we obtained agreement with published work, 
and those cases where we found the published work to be incorrect. In the latter cases, we have isolated 
the error(s) in the published work and report on how we corrected those errors. Section 8 contains a short 
"catalog" of values of 7 for some elementary loop trajectories. This is followed by a short conclusion. 

Note: throughout this paper we use the metric signature (—,+,+,+), and denote Newton's constant 
by G. Prom here on we use units with the speed of light c = 1. 

Section 2: COSMIC STRING MOTION IN FLAT SPACE 

The trajectory of a cosmic string describes a two-dimensional world-sheet in space-time. Points on the 
world-sheet have space-time coordinates given by 

x^=x^'{e,e), (2.1) 

where is a timelike and is a spacelike coordinate on the world-sheet. The string is described by the 
Nambu action, which is proportional to the area of the world-sheet: 

S = -,,j[-g^^^'l^dH. (2.2) 

Here is the mass-per-unit length of the string, 5^^^ is the determinant of the 2-dimensional metric on the 
world-sheet induced from the Minkowski metric, and the integration is over the entire world-sheet of the 
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string. If we define a;**, „ = where a = 0, 1, then the induced 2-diniensional metric is given by 

5i6^=5^-a;^„x^6. (2.3) 

If we denote time and space derivatives on the world-sheet by ' = ^ and ' = ^ , then the determinant g^'^^ 
is 

^ — X X jj X ^ X jj^X X 1/ • ^2.4^ 

Note that x'^ is spacelike and is timelike. 

The Lagrangian equations of motion for the string are rather cumbersome [3]. However, the action 
(2.2) is invariant under the reparameterization (gauge transformation) — > ^''(^), so the equations can be 
simplified by a judicious choice of the parameters One may choose the parameters so that a;'* satisfies 
the gauge conditions 

x^'x'^ = and x^^x^ + x'^x'^ = 0. (2.5) 
With this choice of gauge, the equation of motion is the two-dimensional wave equation 

_ yjiii = 0. (2.6) 

The gauge conditions (2.5) still allow a further reparameterization where ^ = and ^'^ = ^ . Together 
these imply that = This allows us to set ^ = t. If we rename = a, then the coordinates of the 
string world-sheet (2.1) become 

x^" = x^'{t,a). (2.7) 
With this choice of parameters, the gauge conditions (2.5) become 

x^x[ = 0, and x'xi + x'^x\ = 1, (2.8) 

where the index i = 1,2,3 is a spatial index. The equation of motion becomes the two-dimensional wave 
equation 

- x"' = 0. (2.9) 

The time part of the equation of motion (2.6) is satisfied automatically. 
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The general solution to the equation of motion (2.9) is 




(2.10) 



Here, the function a defines the left-moving and b the right-moving component of the string. The first gauge 
condition applied to x implies that a'^ = b'^, where here the prime means differentiation with respect to the 
function's argument. The second gauge condition implies that a'^ + b'^ = 2. Together the gauge conditions 
force the functions a and b to satisfy 



Up to this point, our treatment of cosmic strings includes both the case of infinite strings and the case of 
closed string loops. From here on, to study gravitational radiation, we consider only the case of closed loops. 
In this case, the world-sheet of the (assumed non-self-intersecting) string has the topology of a cylinder 
M X S^, and will be referred to as a "world-tube". Because the string forms a closed loop, one finds an 
additional constraint on the otherwise arbitrary functions a and b. 
If the cosmic string has the form of a closed loop, it follows that 




1. 



(2.11) 



x{t,a + L) = x{t,a) ^ cr,t, 



(2.12) 



where the constant L is the length of the loop. This implies that 



a{t + a) + b{t - a) = a{t + a + L) + b{t - a - L) V cr, t. 



(2.13) 



If we define the null coordinates u and v by 



u = t + a, V = t — a, 



(2.14) 



then (2.13) becomes 



a{u + L) — a{u) = b{v) — b{v — L) M u, v. 



(2.15) 



However, because u and v can be varied independently, it must be the case that 



a{u + L) — a{u) = b{v) — b{v — L) = c, 



(2.16) 
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where c is a constant vector. If we choose to work in the center-of-mass frame of the loop, then c = 0. This 
follows since in the center-of-mass frame we have 

L L 



= j Sdn ^ J ^[a'{t + (j) + b'{t- a)]da 



(2.17) 





= ^[a{t + L)- b{t - L) - a{t) + b{t)] = c. 
Thus, in the center-of-mass frame, the functions a and b are periodic with period L, 

a{t + a + L)=a{t + a), b{t - a - L) = b{t - a). (2.18) 

Because the functions a and b are periodic, each can be described by a closed loop. These loops will be 

referred to respectively as the a-loop and the b-loop. Together, the a- and b-loops define the trajectory of 
the string loop. 

Because the functions a and b are periodic in their arguments, the string loop is periodic in time. The 
period of the loop is L/2 since 

x{t + ^, (7 + ^) = ^[a{t + (7 + L) + b{t- a)] = ^[a{t + a) + b{t - a)] = x{t, a). (2.19) 

For the remainder of this paper, we will set the loop length L = 1. The period of the loop is then At = 
L/2 = 1/2, and the section of the world- tube swept out by the loop in a single oscillation is covered by the 
coordinates a e [0, 1] and t £ [0, 1/2]. The entire world- tube is covered by u e [0, 1] and t € (— oo, oo). 

The reason that one may set L = 1 is remarkable: the power radiated in gravitational radiation from a 
loop of a given shape is invariant under a rescaling (magnification or shrinking) of the loop, provided that 
the velocity at each point on the rescaled loop is unchanged [3]. A formal proof of this is given in [14]. Thus, 
to calculate the radiated power it is sufficient to consider only those loops with total length L = 1. 

The null coordinates u and v defined in (2.14) are more convenient than the coordinates t and a. The 
u, V coordinates are called null because the tangent four- vectors duX" and d^x" associated with them are 
null. This follows because 

d da d dt d d d ^ 
du du da du dt da dt ' 



(2.20) 



d _ da d ^ dt d _ d ^ d 
dv dv da dv dt da dt ' 



are null vectors in Minkowski space. To be more explicit, since points on the world-tube have space-time 
coordinates [t, x] given in terms of u and v by 

x°'{u,v) = ^[u + v, a{u) + b{v)], (2.21) 

the gauge conditions now imply that 

duX^'duXa = ~(^] a'^ = 0, and 

^ \ ^ 2 (2-22) 
d,x"d,Xc, = -(^^^ b'^=0. 

Because of the periodicity of the loops, the world-tube may be covered by the coordinates u and v in many 
equivalent ways. One convenient covering is to take u € [0, 1] and v £ (— oo, oo). The region of the world- 
tube swept out in a single oscillation of the loop is then covered by u e [0, 1] and v G [0, 1]. Note however 
that this is not the same region of the world-tube as t e [0, 1/2] and a e [0, 1]. This is shown in Figure 1 of 
reference [14]. 

The energy-momentum tensor T'*'' for the string loop may be found by varying the action (2.2) with 
respect to the metric. In flat space, with our choice of coordinates and gauge, it is 

1 CXD 

T''^(y«) =fj,J du J dv G^''{u,v)5^{y" - x"{u,v)), (2.23) 

-oo 

where G^"^ is deflned by 

G^''(m, v) = dux^^dyx" + dyx^'dux''. (2.24) 

Note that the volume element for the (m, v) coordinates is related to that of the coordinates {t, a) by the 
Jacobian of the coordinate transformation. Thus, 

dudv = 2dadt. (2.25) 

Because of the delta function which appears in (2.23), the stress-tensor T'^'^ vanishes everywhere except on 
the world-sheet of the string loop. 
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Section 3: POWER RADIATED IN GRAVITATIONAL RADIATION 

The power emitted by an oscillating loop in the form of gravitational radiation may be determined in 
the weak-field limit. This is an excellent approximation for cosmologically interesting cosmic strings because 
the amplitude of the metric perturbation /i^,^ is of order G/x/c^ w 10~^. Because the gravitational radiation 
is weak, its back-reaction on the loops does not modify a loop's motion significantly in a single oscillation. 
Hence we calculate the rate of gravitational radiation in the approximation that the back-reaction can be 
neglected, so that a loop oscillates periodically in time. 

The standard formulae used to calculate the power lost to gravitational radiation typically assume that 
the energy of the source is gradually dissipated into radiation, and that the stress-energy tensor of the source 
vanishes with time. In our case, the source is a periodically oscillating loop whose stress-energy tensor does 
not vanish with time, and therefore the standard formulae require minor modifications. Since the loops that 
we study in this paper have period 1/2, they radiate only at discrete angular frequencies given by: 

LVn = 4:TTn for n = 1, 2, 3, • • • . (3.1) 

The power radiated per unit solid-angle into the nth mode is given by the standard formula (equation 10.4.13 
of [15]) 

HP C 1 

^ = -a;^[T;,(a;„0)T''^(a;„f2) - -\r\{u;r,n)f]. (3.2) 
In this equation, the Fourier transform of the stress-energy tensor is defined by 

Ti^^un) = lim 1- r dt I d'xe'"^'-''-^%4t,x), (3.3) 

T->00 ZI J_rp J 

and * denotes complex conjugation. The vector symbol over x denotes an ordinary flat-space three-vector, 
and O is a unit-length three-vector with spatial Cartesian components given by 

Q= {cos (j) sin 9, sin (j) sin 9, cos 9). (3.4) 

To calculate the total radiated power, one must integrate over all directions on the unit sphere. This 
introduces integrals of the form 

dnf{n)= sin9d9 d<f)f {9 , </)) (3.5) 
Jo Jo 
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into the equations that follow. 

The Fourier transform of the stress-energy tensor (3.3) is defined as the limit of a infinite-time integral. 
This differs slightly from the case of non-periodic sources. It is easy to see that T^,/(a;f2) vanishes unless ui 
takes on one of the discrete values a;„. (This is shown in reference [14] following equation (3.11)). For these 
values of a;„, the T — > oo limit of the integral is not hard to calculate: since the source is periodic, the limit 
appearing in (3.3) is equal to the integral over a single period, i.e., the same integral with 2T set equal to 
1/2. Substituting expression (2.23) for the stress-energy tensor into the formula for t^v one finds that 

T^^{u>nfl) = 2fj f du f d?;G^^(u,'y)e''^"("+''-"-("(")+^(''»)/2_ (g g) 

Jo Jo 

Note that the limits of integration in (u, v) space cover one complete oscillation period of the world-sheet of 
the string loop. 

The total radiated power is obtained by summing over all the modes: 

n=0 

The n = mode has been included in the sum for later convenience; it makes no contribution because of the 
factor of appearing in (3.2). This expression for the total power can be put into a more useful form by 
using the explicit formula (3.6) for the Fourier transform of the stress-tensor. This gives 

P=^^ y ul [ dnC du C dv C du /'d{)^(M,^;,ii,{))e*'^"(^*(™")-"-^^(™^^^ (3.8) 

J Jo Jo Jo Jo 

where we have defined 

V(w, V, u, v) = G^^iu, v)G'"'{u, v) - ^G\{u, v)G\{u, v). (3.9) 

The functions At = {u+v — u—v)/2 and Ax = {a{u) + b{v) — a{u) — b{v))/2 in (3.8) describe the temporal and 
spatial separation of the two points on the string world-sheet with coordinates (u, v) and (u, v) respectively. 
To save space, in some of the formulae that follow, the arguments of At and are not shown; they should 
be implicitly understood. Since each term in the sum over n equals its complex-conjugate, P is explicitly 
real. For this reason, since the n = term does not contribute, we have changed the sum over n to a sum 
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from — oo to oo at the expense of introducing an overall factor of 1/2 into the formula. Prom here on, this 
sum will simply be denoted by 

It is possible to carry out both the sum over n and the integral over solid angle in closed form. To see 
this, consider the integral 

I{t,r) = Y,^n f dne^^^^-^-^. (3.10) 

n ■' 

The integral over solid angle / rfOe*^'^ is easily evaluated, and equals 47rsin(|i])/|i^. Hence one has 

A-rr 

Because of the absolute value signs that appear, some care is required to obtain this last result - one must 
separately consider both possible signs of a;„. The sum over n may now be explicitly carried out. Using the 
standard formula e*"^ = 2'K6p(0) for the periodic delta function on the interval [— tt, tt), one obtains 

47r2 

I{t, r) = -^('5p(47r(i + |r1)) - 5,(47r(t - |r|))). (3.12) 

Noting further that the periodic delta function may be expressed in terms of the ordinary Dirac delta function 
as 

oo 

5p{x)= ^ 5{x + 2'Kk), (3.13) 

fc— — OO 

and noting that \r\ is always positive, one may combine the two delta functions in (3.12) to give 



I{t,r) = ^ E W + k/2 + \i^)-6{t + k/2-\T\)\ 
i\r\ , ^-^ 



OO 



(3.14) 

= 2m <t + k/2)5{{t + k/2f -\ff), 

k=—oo 

where e{x) = 26{x) — 1 is +1 for a; > and — 1 for x < 0. Prom the definition (3.10) of I{t, f) it is clear that 
applying the derivative operator —id/dt brings down an additional factor of a;„. Inserting the time-derivative 
of (3.14) into (3.8) leads to 

P = AGn^ j du ( dv ( du f dv %lj{u,v,u,v)^e{At + k/2)6{{At + k/2f - {Axl"^). (3.15) 

ft=_oo-^o Jo Jo Jo oAt 

In this expression, the dependence of At and Aa; on the four variables u, v, u and v has not been explicitly 
shown. The derivative operator means: first take the derivative of I{t,r) with respect to the first 
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argument, then substitute in the functions At and Ax. Alternatively, it refers to any combination of derivative 
operators in (3.8) which will bring down a factor of iUn- 

It is possible to re-express the sum over k of these four integrals as a single integral, simply by shifting 
one of the integration variables to the range — oo to oo. For example, if we choose to shift v, then because the 
functions a and b are periodic, it is easy to see that Ax{u, v,u,v — k) = Ax{u, v, u, v) and that ip{u, v,u,v — 
k) = tp{u,v,u,v). However, the time function is noi periodic; one has At{u,v,u,v — k) = At{u,v,u,v) + k/2. 
Since the period of the loop is 1/2, the energy radiated in a single oscillation of the loop is thus given by 
E = P/2: 



Note that the choice to shift v was arbitrary; we could have chosen to shift any one of the four integration 
variables. If we had chosen to shift some other variable, the only changes to (3.16) would be that the 
integration range for v would be from to 1, and the integration range of the new shifted variable would be 
from — oo to oo. This expression for the energy radiated into gravitational radiation during one oscillation 
of the cosmic-string loop can be evaluated exactly in the case of piecewise linear string loops. 

To make this formula directly useful, one must replace the operation d/dAt by an explicit operation in 
terms of derivatives w.r.t. the variables u, v, u and v. The desired effect of this operation is to bring down a 
factor of iijJn when applied to exp(ztt'„[At(u, v,u,v) — Q. ■ Ax{u, v, u, v)]). Let us denote this operation by 



E = 




(3.16) 



D{u^ u^v) = U {u^ v)du + V{u, v)dv — U (w, tt, v)du — V{u, v)dy 



(3.17) 



where the functions U, V, U, and V are determined by the desired effect of D on the exponential: 



D exp{iu)n[At — ■ Ax]) = iUn exp{iUn[At — fl ■ Ax\). 



(3.18) 



Because D is chosen to be a linear differential operator, (3.18) is equivalent to the four equations 



DAt{u,v,u,v) = 1 and DAx{u,v,u,v) = 0. 



(3.19) 



Substituting in the definitions of At and Ax, this may be written as the 4x4 matrix equation 



1 1 1 1 H ^ 

d'{u) b'{v) a'iii) b'{v) ) U 

\V 




(3.20) 
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The solution to these linear equations yields the following expression for the differential operator D: 

D{u, V, u, v) = 

J'{v) ■ [a'jn) X b'{i,)]du - a'ju) ■ [a'{u) x b'{v)]d^ - a'ju) ■ [b'jv) x b'iv)]di + a'ju) ■ [b'{v) x a'iu)]d^ 
P{v) ■ [a'{u) X b'{i)] - a'iu) ■ [a'({t) x b'{i)] + a'{u) ■ ^'{v) x b'(v)] - d'{u) ■ x a'({t)] 

(3.21) 

In terms of the operator D, the energy radiated in gravitational waves during one oscillation of the string 
loop is given by 

E = 2Gii^[ du I dv I du I dvi}{u,v,u,v)D{u,v,u,v)(e{At)5{{Atf-\Ax\^)\. (3.22) 

Jo Jo Jo J-oo \ ) 

In this expression the differential operator D acts on all the quantities that stand to its right. 

Before continuing, we note that half the terms may be easily eliminated from (3.22). Under the operation 

of interchanging the variables (m, w) with («,«), U{u,v,u,v) = U{u,v,u,v), V{u,v,u,v) = V{u,v,u,v) and 
At{u,v,u,v) = —At{u,v,u,v), while ■^(u, v) = tjj{u,v,u,v) and the arguments of the 6 function are 
invariant. Using these results, we interchange the variables {u,v) with {u,v) in the first two terms of E 
in (3.22). Recalling that we could have shifted any one of the integration variables in (3.15), as explained 
following (3.16), one finds that 



(3.23) 



/•I /•! rl roo / \ 

du dv du dv^p{u,v,u,v){Udu + Vdy)(e{At)5{{Atf-\Ax\^) 

Jo Jo Jo J-oo V / 

= [ du [ dv [ du [ dvi;{u,v,u,v){Udu + Vdy)(e{-At)6{{Atf -\Ax\^)). 

Jo Jo Jo J-oo V / 

This equation, along with e{—x) = —e{x), may now be used in (3.22) to eliminate two of the four terms in 
the operator D, yielding 

E = 4Gn^[ du f dv f duf dv'4){u,v,u,v){Udu + Vdy)(e{At)5{{Atf-\Ax\'^)\. (3.24) 

Jo Jo Jo J-oo \ J 

Although (3.24) has only half as many terms as (3.22), it is not the most useful form for our purposes. 

The most convenient expression for E is obtained by replacing e{At) in (3.22) by 26{At). The e{At) = 
20{At) — 1 term can be replaced by 29{At) because the —1 term makes no contribution to the integral. 
To see this, consider the efi'ect of replacing e{At) by —1 in (3.22). Denoting this by -E(-i) and again using 
the transformation properties of U{u,v,u,v), V{u,v,u,v), ip{u,v,u,v) and the arguments of the 5 function 
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under the operation of interchanging the variables {u, v) with (u, v), one finds that 

E(_i^=-2Gfi^ [ du f dv f du I dv%i}{u,v,u,v)(udu + Vdy-i}du-Vd^5{{Mf~\Ax\^) 

Jo Jo Jo J-oo \ J 

= -2G/i^/' du [ dv I du I dt5 t^, u, ^) f + - C/as - yac)<5((Af)2 - |Af = 0. 

Jo Jo Jo J-oo \ J 

(3.25) 

This proves that -B(-i) = and thus that e{At) in (3.22) can be replaced by 26{At). 

The integration range v G (— oo, oo) in the expression for E (3.22) may be replaced by the finite range 

V G [—2,2]. This is because the integrand of (3.22) vanishes unless v € [—2,2]. Physically this is because 
in the center-of-mass frame, the string loop remains centered about a fixed coordinate location. Since the 
6 function has its support only on the light-cone, it is impossible for regions of the string located far in the 
past or future to interact. To see this result mathematically, first consider the argument of the 6 function. 
The only contributions to E arise when this argument is positive, i.e., when 

u + v-u-v>0. (3.26) 

Since all three integration variables u,v,u lie in the range [0, 1], the ^-function vanishes unless the variable 

V lies in the range v e (— oo, 2]. Further restrictions arise from considering the argument of the 5 function. 
The only contributions to E arise when this argument vanishes, i.e., when 

{u + v-u-vf = [a{u) - d{u) + b{v) - b{v)f. (3.27) 

However since the total length of the a-loop is 1, the maximum length of the vector a{u) — a{u) is 1/2. 
Similarly, the maximum length of the vector b{v) — b{v) is 1/2. Hence the largest possible value of the r.h.s. 
of (3.27) is 1. This shows that the integrand vanishes unless the quantity 

u + v -u-v e[-l,l]. (3.28) 

Again making use of the possible ranges of u, v, and u this implies that the 6 function vanishes unless v G 
[—2,3]. Combined with the restrictions arising from the 9 function, this implies that the only contributions 
to E arise from the range v G [—2, 2]. Thus we obtain the final form of our result 

E = AGi/( du I dv I du I d{5 ^/'(w, {;)£>(«, u, u, ti) f6»(At)(5((Ai)2 - [Axj^)). (3.29) 

Jo Jo Jo J-2 V / 
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The expression just obtained has an interesting physical interpretation. 

The string loses energy to gravitational waves precisely because of the gravitational self-interaction of 
the string with itself. From this point of view, the integral over {u,v) in (3.29) is an integral over the 
sources (the stress-tensor of the string worldsheet at x°'{u,v)) that create a metric perturbation at the 
space-time point x°'{u,v). The metric perturbation at x°'{u,v) is obtained by multiplying the source times 
a retarded propagator and integrating over the entire history of the source: the part of the world-sheet that 
can contribute is covered by the coordinates u G [0, 1] and v G [—2,2]. The metric perturbations from the 
loop at one space-time point propagate along the light-cone from that point to interact with some other 
point on the loop at some later time. The product 65/2-1: that appears in (3.29) is precisely the retarded 
propagator (equation 12.133 of [16]). This creates the mechanism for energy loss: the string must do work 
against the tidal forces created by the metric perturbations due to the string itself. The energy lost in a 
single oscillation is obtained by integrating this work over the region on the string's world-tube covered by 
the coordinates u G [0,1] and v G [0,1]. Thus, the loss of energy due to gravitational radiation may be 
thought of in terms of a loop which creates metric perturbations, interacting with a loop whose motion does 
work against these perturbations. 

We will see in the following sections that (3.29) can be evaluated exactly in closed analytic form for any 
piecewise linear cosmic string loop. The value of 7 is then given immediately by (l-l)- 

Section 4: PIECEWISE LINEAR LOOPS 

We now restrict our attention to piecewise linear loops. These are loops for which the functions a{u) 
and b{v), which define the loop's trajectory, are piecewise linear functions. The functions a{u) and b{v) may 
be pictured as a pair of closed loops, which consist of joined straight segments. The segments on the a- and 
b- loops join together at kinks where d'{u) and b'{v) are discontinuous. The a- loop has Na linear segments, 
and the b-loop has Nb linear segments. Part of a typical a- loop is shown in Figure 1. 

The following conventions, also shown in Figure 1, are used to describe piecewise linear loops. The 
coordinate u on the a-loop is chosen to take the value zero at one of the kinks, and increases along the loop. 
The kinks on the loop are labeled by the index i where i = 0,l,2,...,A^a — 1. The value of u at the ith kink 
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is denoted by Ui. Without loss of generality we set uq = 0. The segments on the loop are also labeled by the 
index i; the ith segment is the one lying between the iih kink and the {i + l)th kink. The kink at w = u^^ 
is the same as the first kink at u = wq = 0. Since |a'p = 1, the a-loop has length 1. Thus, even though uq 
and ujv,. are at the same position on the loop, uq = while ujv„ = 1- The entire range u e (— oo, oo) may 
be covered by allowing the coordinate u to continue around the a-loop in a periodic way. This also extends 
the range of the index to i € Thus, for example, U-jv„ = — 1, mq = 0, mjv„ = 1, M2JV„ = 2 and so on are all 
located at the same position on the a-loop. Because |a'(u)| = 1, the length of the linear segment between 
the kinks at Ui and Ui+i (the iih segment) is 

Aui = Ui+i - Ui. (4.1) 

The loop's position a{u) at u = Ui is denoted 

cii = a{ui). (4.2) 

The constant unit vector tangent to the ith segment (pointing in the direction of increasing u-parameter) is 
denoted 



Aui 

With these definitions, the function a{u) for u € [ui,Ui+i] may be written 



(4.3) 



d{u) = di + a'i{u - Ui) for u € Wj+i]. (4.4) 
Note that for consistency, putting u = Uj+i in (4.4), one must have 

di+i = di + a[{Aui). (4.5) 
Similar notation is used for the function b{v). For v € [vj,Vj+i], the function b{v) may be written 

b{v) = bj + b'j {v - Vj ) for u e [vj , Vj+i] . (4.6) 

Thus the a-loop is entirely specified by the quantities di; from them one can obtain both Au, = \di+i — di\ 
and a- given by (4.3). Identical notation is used for the 6-loop. 
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The function il>{u, v, u, v) takes a special form in the case of piecewise Unear loops. In general one has 

2 (4.7) 

^^diiX^dyX dy^X jj^dy X 1/ ~\~ dy^X^ QyX dyXi/dyX ^ dyX dyX\dyX^ dyXi^^ 

where, again, x°' = x"{u,v) and = x°'{u,v). For the purpose of evaluating integral (3.29) it is necessary 
to break up the integrations over {u, v,u,v) into rectangular four-cells. Each four-cell is denoted by a set of 
indices {i,j,k,l). The indices refer to segments on the a- and b-loop, each of which defines a specific range 
for one of the coordinates (u, v,u,v). The index i will always refer to the u coordinate, j to the v coordinate, 
k to the u coordinate and I to the v coordinate. Within each cell, we may write the four-vectors and tangent 
vectors 



x" {u,v) = i [u + i;, flj + a ■ (li — Uj) + bj + 6' {v — Vj 



= i[l,a^], a,x" = i[l,6'] 



for u e [ui, Ui+i] and v e bj, Vj+i] (4.8) 



and similarly 



x°'{u,v) = ^[u + v,ak + d'f.{u ~ Uk) + bi + b'i{v - vi) 
dux^ = ^[l,ai], dtx^ = ^[l,b[] 



for it e [uk,Uk+i] audi; e [vi,vi+i]. (4.9) 



Using (4.8) and (4.9) in (4.7) we find that ip{u, v, u, v) is a constant, tjjijki , when (u, v, u, v) are in the intervals 
defined by the segments k, I). For any set of segments k, I), the constant tpijki is given by 

^ijki = lii-l + a[ ■ 4)(-l + b'^ ■ b[) + (-1 + d[ ■ b[){-l + b'^ ■ 4) - (-1 + d', ■ ^.)(-l + b\ ■ 4)]. (4.10) 

Note that ipijki G [—1,5/4]. Also note that ipijki vanishes when {i — k) mod A^a = or (j — I) mod = 0. 

It is helpful to keep track of whether the indices in a given equation refer to kinks or segments on the a- 
and b-loop. For instance, w, refers to the value of the parameter u at the ith kink on the a-loop. Similarly, 
Bi denotes a vector from the origin to the kink at u = Ui on the a-loop. By contrast, a- is a unit vector 
parallel to a specific segment on the a-loop. The indices on Vijfc; also refer to specific segments on the a- 
and b-loops. The segments define specific ranges for the coordinates {u, v, u, v). Of course these ranges will 
change as the indices take on different values. 

The formula (3.29) for the energy radiated in gravitational waves during one oscillation of the string 
loop may now be rewritten for the case of a piecewise linear loop. The integrals over {u, v, u, v) in (3.29) may 
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be broken up into a sum of integrals over the individual segments making up the a- and b-loops. Because 
•ipijki is a constant in each of these integrals, it may be pulled out of the integration, giving 

E = 4Gii^J2 H £ I] ^^i-^' du dv du dvDijki{u,v,u,v)[e{At)d{{Atf -\Ax\^)]. 

'=0 k=0 l=-2N, I J. 

(4.11) 

Note that the summation of I is not from — oo to oo but is only over the finite range corresponding to 
V e [—2,2] as shown following equation (3.28). Here 

Dijki(u,v,u,v) = Uijkidu + Vijkidy - ilijkidu - Vijkidv, (4-12) 

where the coefficients of the derivative operators are constant on any k, I) segment, and are given by 

Uijki = Qijki b'j ■ [a'fe X b'i\ Uijki = Qijki a[ ■ [6^- x h'l] 

(4.13) 

Vijki — —Qijki a- • [a'j. X h'l] Vijki = —Qijki d[ ■ [b'^ x a^], 
with (twice) the inverse determinant given by 

Qm ='^(b'j- [4 X b'l] - a!, ■ [4 X h'l] + 4 • X h'l] - 4 • X 4]^ . (4.14) 

Note that Uijku Vijki, Uijki and Vijki are all constant for a given set {i,j,k,l). Again, note that the indices 
k, I) in equations like (4.12) refer to specific straight segments on the a- and b-loops. They do not refer 
to the components of some tensor. 



Section 5: EVALUATING THE INTEGRALS 

The four partial derivative operators in Dijki{u,v,u,v) in (4.11) may be trivially integrated over u, v, 
u or V. Carrying out these integrations, E takes the form 

iVa-l Wi,-! JVo-l 2Ari,-l 

i=0 j=0 k=0 l=-2Ni (5-1) 

— fj _ \ _ t7 ( q(^) _ qi'") ^l 

'~^tjkl\'JiJ,k+l.l '^■Lj,k,l) ^tjkl\'JiJ,k,l+l '~'iJ,k,U\' 

where the superscript on each 5* denotes the variable which has been integrated out in that term. The 
quantities •Sj^^p '^Ijli' ^ijli ^ijii appearing in (5.1) may all be expressed in terms of a three-dimensional 

integral containing a 5 function: 

S{Ax, Ay, Az,T, s, C{M, N, a, b, c, rf), A, /J2, p3,p4, p5, p6, P7) 

Ax Ay Az 

= J dx J dy J dz0{T + s{x + y-z))S{f3ixz + l32yz + l33xy + f34Z + (35X + l36y + f37)- 

(5.2) 



000 
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Notice that the Umits of integration have been shifted so that the lower Hmit is always zero. The function 
C is defined by 



C{M, N, a, b, c,d) = < 



' 1 


if 


a—c—l 
M 


_ d-b _ 

N 


integer 


2 


if 


a—c—l 
M 


d-b-1 
N 


= integer 


3 


if 


a— c+1 
M 


d-b-1 
" N 


= integer 


4 


if 


a— c+1 
M 


d-b 
~ N ~ 


integer 



(5.3) 



otherwise. 

The function S does not depend upon C, since C does not appear on the r.h.s. of (5.2). However C provides 
a convenient means to later simplify certain special cases that arise. The quantities S^^j^i, S^^li, S^^jli and 
^ijki given in terms of S by 

5g = S{AviAuk, Avj, Mijki,-1, C{N,, N,, I, k,j, i), {1-b'j- 6j), (1 - 4 • b'j), (4 • b[ - 1), 

{b'j ■ Nijki - M.jki), {M.jki - b'l ■ N.^ki), {M,jM - a'k ' \{N-jki " M^^i)) 



S\% = S{AuuAvuAui, Mijki,-l, C(iV„ iVb, k, I, i,j), (1 - a; • 4), (1 - • b'l), (^ • b[ - 1), 
{3fi ■ Nijki - Mijki), {Mijki - ctk ■ Nijki), {Mijki - b'l ■ Nijki), l{N^jki - Mfjki)) 



Sl% = S{AvjAui, Avi, Mijki,+1, C{N,, N,,j, i, I, k), {1-b'j- b\), {l-a^ b[), (a^ • b'j - 1), 

{Mijki - b'l ■ Nijki), {b'j ■ Nijki - Mijki), {a'i ■ Nijki - Mijki), ^(iV^.^^ - M^;,,)) 



(5.4) 



Sl;li = S{Aui,Avj,Auk, Mijki, +1, C{Na, Ni,, k, I), {I -a',- 4), (1 - 4 • b'^), (a^ • b'^ - 1), 

{Mijki - a'k ■ Nijki), (4 • Nijki - Mijki), {b'j ■ Nijki - Mijki), ^{N^jki - Mfjki)), 
where we have defined 

Mijki =Ui + Vj -Uk-vi, and 

Nijki = ai + bj -CLk-bi- 
Note that the three remaining integrations in S may be done in any order. The relationships 

S{Ax, Ay, Az, r, s, pi,p2, Ps, Pi, /?5, p6, Pr) 
=S{Ay, Ax, Az, t, s, p2,Pi,p3, Pa, Pq, P5, Pi) 
=S{Ax, -Az, -Ay, r, s, -ps, P2, -pi, -Pe, p5, -Pi, Pi) 



(5.5) 



(5.6) 
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allow one to rewrite the integrations in any order. 

It should be noted that in the definition (5.2) of S, there are no terms in the argument of the 5 function 
that are quadratic in x,y or z. This is because the terms u'^,v^,u'^ and in At"^ — \Ax\'^ appear with 
respective coefficients 

-^(l-^-^), -^(1-4-4). -li^-b'rb'i). (5.7) 

It is easy to see that these coefficients all vanish since |a'p = = 1. This is because, as described in 
section 2, the coordinates u and v are null coordinates. 



Section 6: EVALUATION OF S 

The previous section reduced the problem of determining 7 in the piecewise linear case to evaluating 
a set of integrals defined by the function S in equation (5.2). In this section we carry out that evaluation 
in closed form. The 6 function appearing in (5.2) allows us to reduce the number of integrations in 5 from 
three to two. This S function may be written as 5{f{x,y,z)), where the argument of the 6 function is 

f{x, y, z) = pixz + p2yz + (33xy + piZ + P5X + Pey + /J7. (6.1) 

The 6 function will only have support when f{x,y,z) = 0. Solving f{x,y,z) = for x{y,z), y{x,z) and 
z{x,y) respectively, we find that 

_P2yz + P4Z^H3ey±lh , . 

""^y'''- f3,z + (3sy + p, ' ^^-'^ 

Pixz + l3-,x + 34Z + (3r 

PSX + P2Z + P6 

and 

^'^^^ Pix + P2y + P, ■ ^ ' 

The surface z{x, y) consists of a pair of disconnected hyperbolic sheets as shown in Figure 2. The sheets 
are separated by the plane (iix + /32y + /34 = 0, where the denominator of (6.4) vanishes. If these sheets 
pass through the region of integration of (5.2), which is a rectangular box with opposite corners (0, 0, 0) and 
(Ax, Ay, Az), then the 5 function will have support in that region and S may contribute to E. 
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The hyperbolic sheets have a simple physical interpretation. For purposes of clarity, we discuss the case 
^ijki'' other three integrals in (5.4) have similar interpretations. The z coordinate in S parameterizes 
the world-line followed by the fcth kink on the string loop. This kink moves along a straight, null world-line. 
The X, y integrations are over a (diamond-shaped) planar patch of the world-tube, swept out by a linear 
segment of the string loop. Note that these patches are always time-like. The edges of this planar patch are 
bounded by the straight, null world-lines of ith and jth kinks. The future and past light-cones of any point 
on the fcth kink's world line (i.e., fixed z) may intersect the planar patch bounded by the ith and jth kinks. 
Even if they fail to intersect this patch, they will intersect the infinite 2-plane passing through the patch, 
which is parameterized by x and y. The intersection will trace out a hyperbola on the plane which may or 
may not intersect the actual x, y integration region. The hyperbola is given (with z fixed) by (6.3). This 
hyperbola corresponds to the intersection of a 2; = constant plane with the hyperbolic sheets of Figure 2. 
Because one branch of this hyperbola lies on the future light-cone and the other lies on the past light-cone, 
we refer to these as the future branch and the past branch. The two branches are disconnected except in 
the case where the plane passes through the origin of the light cone. In this case they touch at a single 
point - the apex of the light-cone. As the apex moves along the world-line of the kink (i.e., z increases) the 
light-cone sweeps out a region on the x-y plane. This can be seen by taking successive z = constant cross 
sections of the hyperbolic sheets, where the constant ranges from to Az. The region swept out in the x-y 
plane will be bounded by the hyperbolae y{x, 0) and y{x, Az). Note that if we restrict attention to just the 
future branches or just the past branches, then the y{x, Az) hyperbola always lies above (and to the right of) 
the y{x, 0) hyperbola in the x-y plane. Therefore, we refer to these as the "Top" and "Bottom" hyperbolae, 
as shown in Figure 3. The hyperbolae shown in Figure 3 will intersect the x, y integration region only if the 
hyperbolic sheets pass inside the integration box shown in Figure 2. 

There are two useful sets of conditions which can be checked immediately to see if S vanishes. These test 
whether the hyperbolic sheets z{x,y) pass through the integration box. If the sheets do not pass through the 
box, then the 5 function in (5.2) has no support, and S vanishes. One of the conditions applies for s = +1. 
In this case, the 9 function in (5.2) restricts the integration to be over only the hyperbolic sheet swept out by 
the future light-cone. The other condition applies for s = —1. In this case, the function in (5.2) restricts 
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the integration to be over only the hyperbolic sheet swept out by the past light-cone. If s = +1, then 

!T + Ax + Ay <0,ov 
f33AxAy + P5Ax + (3eAy + l3T>0,or (6.5) 
Az<T and l3iAz + /Jy < 0. 

Similarly, if s = — 1, then 

!r + Aa; < 0, or 
/34Az + /?7 > 0,or (6.6) 
Ax + Ay<T and PaAxAy + /JgAx + PeAy + /Jt < 0. 
In practice, these conditions are frequently satisfied, so their implementation saves large amounts of com- 
puting time. 

It is straightforward to do the integral over z in (5.2) to eliminate the 6 function in S. We make use of 
the standard formula 



/ 



Sif{z))g{z)dz = ^g{z,)/\^{z,)\ (6.7) 
p 

where the sum is taken over all the roots Zp of f{z) that lie within the range of integration. Because 

^ = Pix + p2y + 04, (6.8) 

the integration of (5.2) over z yields 

S{Ax,Ay,Az,T,s,C,0u---,p7) 

Ax Ay 

_ [ , [ , ejzjx, y))9{Az - z{x, y))9{T + s{x + y- z{x, y))) (6.9) 



The first two 9 functions in this equation arise from (6.7): the only roots included in the sum over p are 
those lying within the range z G [0, Az] of ^-integration. 

The integral (6.9) has a simple physical interpretation which is directly connected to the physical 
interpretation of (5.2) already given following equation (6.4). The x and y integrations are over a planar 
patch on the string loop's world-tube, as in (5.2). The first two 9 functions in (6.9) only have support 
between the hyperbolic curves y{x,0) and y{x,Az) (i.e., the Bottom and Top curves) shown in Figure 3. 
The curve y{x, 0) describes the intersection of the hyperbolic sheet shown in Figure 2 with the bottom of 
the integration box. The curve y{x, Az) describes the intersection of the hyperbolic sheet with the top of 
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the integration box. Thus, the first two 6 functions in (6.9) will only have support if the hyperbolic sheets 
pass through the box of integration in (5.2). The third 6 function in (6.9) effectively restricts the integration 
to be over the future branches (s = +1) or the past branches (s = — 1) of the hyperbolae. 

When evaluating (6.9), there are five different fundamental types of integrals that can arise depending 
upon the relative positions of the two hyperbolic curves y{x, 0) and y{x, Az) and the rectangular region in 
the x-y plane bounded by the opposite corners (0,0) and (Ax, Ay). Each of the five possibilities are shown 
(using the future branches of the hyperbolae) in Figure 4. Each type of integral may also occur with the past 
branches of the hyperbolae. The first type occurs when both hyperbolae (Top and Bottom) pass through 
the planar patch. The second type occurs when the entire planar patch lies between the two hyperbolae. 
A third possible type occurs when neither hyperbola passes though the planar patch and the planar patch 
does not lie between the two hyperbolae. For this type of integral, (6.9) has no support and vanishes. The 
last two types arise when one of the hyperbolae (Top or Bottom) passes through the planar patch but the 
other does not. 

Remarkably, each of the fundamental integral types can be done analytically in closed form. To assist in 
this process, the x integration in (6.9) must be divided into consecutive ranges; the type of integral in each 
range is different. A systematic method for determining these integration ranges in x will be given next. 

The dividing points (in x) between the successive ranges are determined by the four points at which 
the bottom and top hyperbolae y{x, 0) and y{x, Az) intersect the lines y = and y = Ay. We denote the x 
coordinates of the four intersection points by 0i, . . . , The type of integral being done in (6.9) will change 
at each of the intersection points which is within the boundaries of the x integration, i.e. < a; < Ax. The 
boundaries of the x integration will be labeled 0o and 05- The x coordinates 0i, . . . , (/14 are given by; 
.^i=x(0,A^) = -^4^±^ if C 5^1 else </.i = 0, 

PlAz + P5 

(/)2 =x{Ay, Az) = A _L A _L « ^ ^ 2 else (j)2 = 0, 

(3,Ay + p,Az + P, ^^^^^ 

cl>s =xiAy,0) = YJ if C^3 else .^3 = 0, 
PsAy + 

07 

=x(0, 0) = if C 7^ 4 else (^4 = 0. 

The values of C are checked because there are four special cases. In these special cases the formula given 
for one of <;6i, . . . , 04 is indeterminate because the numerator and denominator in (6.10) vanish. 
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The conditions for the four special cases are given in the definition (5.3) of C. Physically, these special 
cases arise when one end of the kink's world-line, parameterized by z, touches one of the four corners of the 
diamond-shaped patch of the world-sheet defined by the x and y integrations. The intersections of the future 
and past light cones of a point on the kink's world-line with the plane defined by the diamond-shaped patch 
are usually hyperbolae. However, when the point on the kink's world-line is also a corner of the diamond- 
shaped region, then the plane passes through the apex of the light cone and the hyperbola degenerates into a 
pair of straight lines. These straight lines will lie along the two edges of the planar patch which are joined at 
the corner where the kink's world-line touches. For each of the four special cases, one of the formula for (j) in 
(6.10) would become indeterminate. Consider, for instance, the case where the end of the kink's world-line 
(ie. z = Az) touches the lower left corner of the integration region (ie. x = y = 0). At this point, we have 
x^lujv) = x"{u,v). In this case the future top curve lies along the left and lower sides of the rectangular 
integration region. One can verify that in this case, both the numerator and denominator in the equation for 
(f)! vanish. It is (pi that becomes indeterminate in this case because (f)i is the x coordinate of the intersection 
of the top curve and the line y = 0, which does not have a unique solution in this case. The other three 
special cases are similar. In each case a different </> would become indeterminate if the value of C were not 
checked. 

All four special cases are dealt with in the same manner. The purpose of the ^'s is to locate the x 
coordinates where the type of integral being done changes. However, since intersection curves which lie 
along the edges of the integration region never cause the type of integration being done to change, it is 
sufficient to simply set the corresponding (f)i to zero (or to any value outside the range {(j)o, (ps))- 

The support of the x integration in (6.9) may be less than the range < a; < Aa; because of the third 
6 function. Thus one does not always have (po = and ^5 = Ax. Because the integrations in (6.9) have 
support only between the future branches (s = -|-1) or the past branches (s = — 1) of the hyperbolae, it is 
convenient to define ^0 and ^5 in a more general way. First, we define to be the vertical asymptote to 
the bottom hyperbola when s = +1 and zero otherwise, and (pr to be the vertical asymptote to the top 
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hyperbola when s = — 1 and zero otherwise, 



a ^ o (6-11) 
=e{-s)x{y oo,Az) = -0{-s) . 

The boundaries of the x integration are then defined to be 

4>o =max(0, 0b) 

(6.12) 

(f)5 = min(Aa;, (pr + 9{s)Ax). 
This definition of (f)o and ^5 eliminates regions of integration which only contain the past branches of the 

hyperbolae when s = +1 and regions which only contain the future branches of the hyperbolae when s ~ —I. 

One may now express 5 as a sum of integrals over the successive ranges of the x-integration. The set 

of points {(1)0 . . .(p5} partition the a;-integration into at most five ranges. Let {a;o, x\,X2-iX:i,Xi, x^} be the 

increasing sorted set of <j)'s, 

{xo,Xi,X2,X3,Xi,X5} = SOrt({0o,01,'A2,</'3,04,05}), (6.13) 

sorted so that Xi < Xi+i. If we define the mid-point between two successive a;'s as 

Xn = {Xn + X„+-i)/2, (6.14) 

then we may rewrite the S integral (6.9) as a sum of integrals over each successive x range; 

4 

S'(Aa;,A2/,A2,r,s,C,/3i,---,/37) = ^9(xn - (/>o)6'(</>5 - 5n)I'(a;n, a;„+i, Ay, Az, r, s, /3i, • • • , /Jt). (6.15) 

The two 9 functions ensure that only the ranges of x between and 05 may contribute. The function T is 
defined to be 

Ay 

rrt A A R R\ [a ^(^(a;, y))9{Az - z{x, y))9{T + s{x + y - z{x, y))) 
T{xi,Xu,Ay,Az,T,s,Pi,---,P7) = J ^x J dy _^ ^^j —■ (6.16) 

xi 

The reason that we have defined the function T is that the x integration being done in (6.16) is over a region 
which contains only one of the five possible fundamental types of integrals discussed in the paragraphs 
preceding equation (6.10) and shown in Figure 4. 

What remains is to find the analytic form of (6.16) for each of the five possible types of integrals that 
can arise. Recall that the different types of integrals arise, as shown in Figure 4, from the different possible 
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relative positions of the Top and Bottom hyperbolae compared to the {x, y) region of integration. If the 
region of integration is between the two hyperbolae, then the first two 6 functions in (6.16) will have support 
over the entire region, and the limits of the y integration will run from to Ay. If both hyperbolae pass 
through the region of integration, then the first two 9 functions in (6.16) restrict the y limits of integration to 
run from y{x, 0) to y{x, Az). The last two types of integrals which give non-zero contributions are when one 
of the hyperbolae passes through the region of integration but the other does not. These types of integrals 
will have y limits of integration which run from to y{x, Az) or from y{x, 0) to Ay depending on whether it 
is the Top or the Bottom hyperbola that passes through the region of integration. For each type of integral, 
the third function in (6.16) restricts the integration region to be between either the future (s = +1) or the 
past (s = — 1) branches of the hyperbolae. If we make the definitions, 



X ={xi +a;„)/2, 
yb =sy{x,e{-s)Az), 
yt =sy{x,9{s)Az), 

(6.17) 

yi=- e{-s)Ay, 

Vu =d{s)Ay, 

6 =e{pix + sp2yb + Pi), 
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then the integrals in (6.16) can be carried out for each of the five different cases. One obtains 



T{xi, Xu, Ay, Az, r, s, /3i, • • • , /J7) 

L, '^'^ Jo \fiix + P2y + 134 



— { L{xi, Xu, /34 + /32A?/) - L{xi, Xu, Pi, (3i) 



for yt<yb< yi, or 
for yb < yi and t/„ < yt 



l-e(s)v(x,Az)+e{~s)Ay 

dx I dy 

XI J 0{-s)y{x.O) 

sS 



1 



\f3ix + p2y + f34\ 



— Q{xi, Xu, P1P3, PlPa + pspi - P2P5, Pipe - P2P7) 

-L{xi, Xu, Ps, P6 + P20{s)Az) - L{xi, Xu, Pi, P4 + P20{-s)Ay)^ 



for yb<yi <yt< yu 



1 



rvixAz) 

dx I dy . . 

XI Jy(x.n) \PlX + P2y + P4\ 

6 



P2 



L{xi, Xu, Ps, Pe) - L{xi, Xu, Ps, Pe + P2AZ) 



for yi < yb < yt < yu 



^0(-s)y(xAz)+e(s)Ay 

dx I dy 

xi Je{s)y{xfi) 

s5 



\Pix + p2y + P4\ 



— L{xi, Xu, p3, p6 + p20{-s)Az) + L{xi, Xu, Pi, Pa + p20{s)Ay) 

P2 \ 

-Q{xi, Xu, P1P3, Pipe + P3P4 - P2P5, PiPe, - P2P7) 



for yi < yb < yu < yt, or 
for yi <yb< yu and yt < yb 



for yu < yb, or 
for yb<yt< yi- 



(6.18) 



Here, the functions L and Q are the "hnear" and "quadratic" integrals defined by 



L[xi,X2,ci,co) = j dx \n\cix + co\ 

Xl 

(a;+ f^)ln|cia; + co| -x 
. {x2 - a;i)ln|co| 



for CiT^Q 



for ci = 



(6.19) 
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and 



Q{xi,X2,C2,ci,co) = j rfa; ln|c2a;^ + cia; + Col 

xi 

L[xi,X2,Ci,Co) 



= < 



L(xi,a;2,l, ) +L[xi,X2,l, ) + 



Cl+^ c\-4:C0C2 ^ 

2 

(x2 -a;i)ln|c2| 



for C2 = 



for < cf — 4coC2 



1 



+ 



. ^(^l + i;'^2 + ^,0,C2) 



for cf - 4coC2 < 0. 



(6.20) 



Thus T, and hence S has been evaluated analytically for all possible cases. Using the results of this section, 
one can carry out the summations in (5.1) to arrive at a final value for 7; the power radiated in gravitational 
waves by a string loop. 



Section 7: TESTING THE FORMULA FOR 7 AGAINST PREVIOUS RESULTS 

In this section we compare the 7 values given by our formulae to previously published values for a 
large number of loop trajectories. The formulae obtained in sections 5 and 6 were directly implemented by 
computer code. In some cases we find disagreement between our results and those previously published. 
There are, in fact, several cases where conflicting results have been published for the same loop trajectories. 
In the cases where a disagreement was found, we have identifled the errors made in the published work which 
led to the incorrect results. In these cases our formulae give the correct values of 7. We are confident that 
they are correct because, in every case, we have shown our results to be consistent with those given by other 
independent methods. The other methods used to confirm our results were the FFT method of Allen and 
Shellard [14] and/or a corrected implementation of the numerical method used by the original author(s). 

While our formulae handle piecewise linear loops exactly, most of the previous work in this area has 
considered smooth cosmic string trajectories (typically providing analytic expressions for the a- and b-loop). 
To compare the results of our formulae to the published values of 7 for these smooth loops, we calculate 7 
for piecewise linear loop trajectories of approximately the same shape. If the number of segments used {Na 
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and Nb) is reasonably large, then the piecewise linear loop trajectory and the smooth loop trajectory will 
be very similar in shape, and we expect that the values of 7 for the two trajectories will be very close. An 
example of how we generate specific piecewise linear trajectories is given later in this section. The rate at 
which 7 converges as the number of linear segments is increased is also discussed. 

Prior to this work, the only fully analytic closed form solution for 7 for any string loop trajectory was 
given by Garfinkle and Vachaspati [11]. They considered the piecewise linear loops defined by a- and b- loops 
which consisted of just two linear segments each; i.e. Na = = 2. This defines a family of loop trajectories 
which depend on a single parameter (f), the angle between the a- and b-loops. As a function of (I), 7 is given 
by 

■?9 r 9 9 1 

(7.1) 



(A^ 32 



sm 



(l+cos</.)ln(— ^— -) +(l-cos<^)ln(- '^—) 

n + COSffl^ H — COSffl^ 



When calculating 7 with our formulae, we used a- and b-loops with three segments (where the length of the 
third segment was much smaller than the other two). This was necessary in order to prevent a singularity; 
exactly parallel segments cause the determinant in (4.14) to vanish. Equation (7.1) is plotted in Figure 5 
(solid line) along with the 7 values (dots) given by our code for a number of loops with different values of (p. 
Since these are piecewise linear loops, we expect our results to be highly accurate. Indeed, the points plotted 
in Figure 5 showing our results had to be enlarged in order to distinguish them from the plot of (7.1). Thus 
our method completely confirms the results of Garfinkle and Vachaspati. 

The next set of loop trajectories which we use to test our formulae is a three-parameter family of 
trajectories first examined by Burden [10]. The three parameters are L, M and </>, where L and M are 
positive integers and </> is an angle in the range [0,7r]. The Burden trajectories are defined by the a- and 
b-loops: 

^(w) = ~^ — [cos{2'kLu)z + sm{2iTLu){cos(f)X + sini/iy)], 

1 (7-2) 
b{v) = — — [cos(27rMt;)i - sin(27rMt;)x]. 

The b-loop winds M times around a circle in the x-z plane. The a-loop winds L times around a circle whose 

plane is at an angle (p with respect to the x-z plane. The Burden string loops are non-intersecting cuspy 

loops in the case M = 1, L > 1 and (j) not equal to or tt. Burden calculated values of 7 for loops with 

M = 1 and L = 1,2, 3, 5, 15 for several values of ^. The values of 7 for loops with M = 1 and L = 3,5 
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were also calculated by Quashnock and Spergel [13] . Using our formula, we calculated the values of 7 for a 
large number of loops, each of which is a piecewise linear approximation to a Burden loop. In addition, we 
calculated a number of 7 values using the FFT method of Allen and Shellard [14]. Our results for loops with 
M = 1 and L = 3,5 are shown in Figure 6 along with the results of Burden, Quashnock and Spergel, and 
the FFT method. We find excellent agreement among all four sets of results. This also shows that piecewise 
linear loops with fairly small numbers of segments {Na = 16L and Ni, = 16M = 16) can provide excellent 
approximations to smooth loop trajectories and provides further evidence that our formulae are correct. 

Values of 7 for the Burden loops with L = M = 1 have been published by Burden [10] , Vachaspati and 
Vilenkin [9], and Durrer [8]. These results, along with the results of the FFT method and our new method 
are shown in Figure 7. There is excellent agreement between four of the sets of results. However, Durrer's 
results for these trajectories do not agree well with the others. 

To understand why Durrer's results do not agree with the others, we recalculated 7 for these trajectories 
using the same numerical method used by Burden, Vachaspati and Vilenkin, and Durrer. This method 
requires one to calculate the average power radiated by a string loop using the formula previously given in 



where P„ is the average power radiated at frequency a;„ = 4?™ and the integration is over the 2-sphere. The 
details of the calculation of P„ can be found in references [8,9,10]. (Note however the following typographical 
errors in reference [8]. The term J;+i(— /sin0) in (A. 6) should be Ji+i{—lii\n9) and y should be replaced by 
—y in (A. 12) and (A. 13).) Because of the infinite sum appearing in (7.3), one must stop calculating the P„ 
numerically at some value of n, and then estimate the contribution to the sum from larger values of n. Since 
the sum may be slowly convergent, this "tail" may give a significant contribution even when the individual 
Pn are very small. For the L = M = 1 Burden loops with ^ ^ or tt, the tail can estimated with good 
accuracy because the P„ fall off as a power law n~^/^ for large n. Durrer's results are incorrect precisely 
because the contribution of the tail was not included at all. In Figure 8 we show Durrer's original results 
and our calculation of the sum of the first 50 terms of (7.3). Note the agreement between these values. We 
also show the results of the first 50 terms plus an estimate of the tail of the sum along with the results from 



(3.7) 




(7.3) 
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our code. It is clear that when the tail is included, Durrer's results then agree with the results found by 
all other investigators. Thus, again we find that our method is in agreement with the results of previous 
authors. 

Before continuing to compare the results of our new formulae to those in the published literature, we 
take a moment to discuss how piecewise linear loop trajectories are constructed to approximate smooth loop 
trajectories. We illustrate the procedure by explaining how the piecewise a- and b-loops were constructed to 
approximate the L = M = 1 Burden loops considered above. The piecewise linear a-loop was constructed 
by dividing the interval u G [0, 1] into 16 equal segments. This defines 16 values of the coordinate u. These 
values were then perturbed by small random amounts so that pairs of segments on the a-loop would not end 
up exactly parallel. (This is necessary to prevent the determinant in (4.14) from vanishing.) The perturbed 
values of u were then used in the first equation of (7.2) to yield the coordinates of the Na = 16 kinks which 
define the a-loop. Finally, the entire a-loop was translated in 3-space so that the kink with parameter u = 
was positioned at the origin. The piecewise linear b-loops were constructed in a similar manner. Each b-loop 
was constructed to have Ni, = 16 linear segments. In all cases, our values of 7 were within 8 percent of 
previously calculated results with an average diff'erence of less than 3.5 percent. Thus, for the purposes of 
calculating 7, a single wind around a smooth circular path is approximated extremely well by a set of only 
16 linear segments. 

We now examine how the 7 values found for the piecewise linear approximation to smooth loop tra- 
jectories depends on the number of segments used. To determine this dependence, we constructed several 
piecewise linear approximations to the L = M = 1 Burden loops using the procedure given above, each 
with different values of Na and Ni,. The results of four such tests are shown in Figure 9. One can see that 
the values of 7 converge quickly as Na and A^f, increase. It is only for loops with values of (p near and 
180 degrees (where 7 diverges) that a large number of segments are needed to obtain good accuracy. The 
relative errors in four approximations compared to the most accurate approximation (curve D in Figure 9) 
are given in Figure 10. The errors decrease rapidly as the number of segments increases. These errors are 
small and are mainly due to the loops which have 4> close to or 180 degrees. We obtained similar results 
for the L = 3, M = 1 and L = 5,M = 1 Burden loops. Further discussion of how the accuracy of 7 in the 
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piecewise linear approximation of a smooth loop depends on the number of segments Na,Nf, is postponed 
until the end of this section. 

We now continue to compare the results of our formulae to those in the published literature. The next 
set of loop trajectories with which we compare our results is a two-parameter family of loops first studied 
by Vachaspati and Vilenkin [9]. The a- and b- loops which define these trajectories are given by 

a{u) =^[sin(27ru)x — cos(27ru)(cos + sin (j)z)] 

ZTT 

la 

b{v) = — [(— sin(67r?;) — (1 — a) sin(27rt;))x 

27r 3 (74) 

— (— cos(67rt;) + (1 — a) cos(27rt;))y 
- (a(l - a))^/^ sm{ATrv)z] 

where a and are constant parameters, < a < 1 and — tt < (j) < tt. Note that when a ~ 0, these 
trajectories are equivalent to the L = M = 1 Burden loops studied above. These loops have also been 
studied by Durrer [8]. The results found by Vachaspati and Vilenkin, and Durrer are shown in Figure 11 
along with the results of the FFT method and the results of our new method for the case a. = 0.5. We find 
that only the FFT method and our new method are in good agreement. The 7 values given by Vachaspati 
and Vilenkin, and Durrer appear to be too small. In fact, their results are lower than the sum of the first 
300 P„ found by the FFT method (see Figure 11). We take the sum of the first 300 P„ to be a lower bound 
for 7 since continuing the sum to larger n or adding an estimate of the tail (or both) will only increase the 
value found for 7. 

There are two possible explanations for the incorrect results given by Vachaspati and Vilenkin. The 
first possibility is that they incorrectly estimated the tail contribution to the sum in (7.3). Vachaspati and 
Vilenkin claim that the sum in (7.3) is rapidly convergent, with P„ oc for large n. However, we have 
found that the sum is actually much less convergent. For example, the power spectrum for the trajectory 
with a = 0.5 and (j) = 7r/2 is shown in Figure 12. In this case, P„ oc n~^'^^ for n in the range 100 < n < 300. 
By overestimating the convergence of the sum in (7.3), one seriously underestimates the contribution due 
to the tail of the sum. The other possible explanation is that the results reported in [9] are actually for a 
different set of loops than those defined by (7.4). We consider this a possibility not only because the reported 
convergence of the sum (7.3) does not agree with our findings, but also because Vachaspati and Vilenkin 
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include a drawing (Figure 4 of reference [9]) of the loop's shape at two different times during its oscillation. 
However, these shapes do not agree with the shapes given by (7.4). We have confirmed that the loops shown 
in Figure 4 of reference [9] are not the same as the loops given by (7.4), however we have been unable to 
resolve whether the values of 7 reported in reference [9] correspond to the loops defined by (7.4) or to those 
shown in the figure [17]. 

We have calculated values of 7 for the loop trajectories (7.4) using our new formulae for several other 
values of the parameter a. When a = 0, the loop trajectory is equivalent to the L = M = 1 Burden loops. 
Thus, for small a, the loop trajectories (7.4) should be similar to those given by (7.2). In Figure 13 we show 
our results for a = 0.01 (solid line). This is compared to results using the traditional numerical method 
(crosses) and the results for a = (dashed line). The results of our new method agree well with those of the 
traditional numerical method. 

The final string loop trajectories with which we compare our formulae were first given by Garfinkle and 
Vachaspati (equation (2.9) of reference [11]). The a- and b-loops for these trajectories are composed of two 
smooth circular arcs joined by a pair of straight segments. The analytic expressions for the a- and b-loops 
are 

d{u) = [sin((5(u) + 2'Kqu)x — coa{6{u) + 2Tr qu) {cos (j)y + sin (j)z)] 

1 (7.5) 

^(^) =7^ — [sin(/3(v) — 2-Kpv)x — cos{l3{v) — 2Trpv)y] . 
2ttp 

Here, p and q are constants in the range [0, 1], (p is the angle between the two loops, and (3 and 6 are defined 

by 

I3{v) = (1 - p)Tr[-2v] and S{u) = (1 - q)(| + 7r[2u]). (7.6) 

In (7.6), [x] is the greatest integer less than or equal to x. Our results for trajectories with {p,q) =(0.6,0.4), 
(0.4,0.8) and (0.9,0.9) are shown in Figure 14. The results of the FFT method for several trajectories with 
(P)9) =(0.9,0.9) are also shown in Figure 14. Again, we find good agreement between the two methods. 
Garfinkle and Vachaspati do not give specific values but claim that the trajectories given in (7.6) have 7 
values around 100. This is consistent with our results. Durrer [8] has also given values of 7 for some of these 
trajectories. For the three trajectories with parameters {p,q) =(0.6,0.4), (0.4,0.8) and (0.9,0.9) and with 
(f) = 7r/2, Durrer reports 7 values of 19, 26 and 42 respectively. However, because of the errors (explained 
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above) in other numerical results presented by Durrer, we do not have confidence in these values of 7. 
Durrer's results appear to be too low, which would be consistent with leaving off the contribution from 
the tail of the sum in (7.3). The agreement between the FFT method and our new method again gives us 
confidence that our formulae are correct. 

We now return to the question of how accurately the 7 values from piecewise linear loop trajectories 
approximate the 7 values from smooth loop trajectories. In particular, we would like to know how the 
difference between the 7 value of a piecewise linear loop with N = Na + Ni, total segments (7jv) and the 7 
value of the smooth loop that it approximates (70c) falls off as a function of N. Unfortunately, we do not 
know how the difference A(A'') = |7oo — 7Ar| depends on A'' in the general case. However, numerical estimates 
may be made for individual loops with the hope that the results will hold in general. In addition, there is 
at least one case which has been investigated where simple analytic formulas exist for both 700 and 7Ar. 

A detailed numerical investigation of A(A'^) has been carried out for two Burden loop trajectories (7.2) 
with L = M = 1. The numerical values of 7jv have been computed over a wide range in A'' for both loops. 
The value of (p was arbitrarily chosen to be 39° for the first loop, and 111° for the second loop (see Figure 
15). We find that for large A'', both sets of results are well approximated by functions of the foy:mA + BN~^, 
where A and B are constants that depend only upon ^. For the first set of results = 39°), we find that 
7jv ~ 52.01 + 181. 64A^^^ for 60 < A'' < 256. By taking 700 = 52.01, we can find a numerical estimate of 
A(A'^). A similar analysis has been done for the second set of loops, where 7jv ~ 64.49 + 97 A3N~^. Figure 
16 shows a Log-Log plot of A(A'') for both sets of loops. By examining the slopes of the curves in Figure 16, 
we find that in both cases, A(A'') falls off like A''"^ for large N. 

In addition to the numerical investigations of A(A'') given above, there is one case where A(A'') has been 
calculated analytically. In a recent paper, Allen, Casper, and Ottewill have found a simple analytic formula 
for the 7 values of string loops in a particular class [18]. String loops in this class have a-loops which lie 
along a line, and b-loops which lie in the plane orthogonal to that line. In particular, when the b-loop takes 
the shape of an Af(,-sided regular polygon, 7jV(, is given by 




(7.7) 
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When the b-loop takes the shape of a perfect circle, 7 is found to be 

1 — cosa; 

7 = 16 / dx « 39.002454. (7.8) 

Jo X 

Since an A^^-sided polygon becomes a perfect circle in the limit as Nh goes to infinity, the difference between 
equations (7.7) and (7.8) give A{Nb). Figure 7 of reference [18] shows a Log-Log plot of A{Ni,). From this 
plot one finds that, in this case, A{Nb) falls off as A^"^"^. The reason that the 7 values from the piecewise 
linear loops converge to the smooth loop limit faster in this case than in the previous cases is most likely 
because the a-loop in this case is already piecewise linear. While we do not know how the errors in 7jv scale 
with increasing A'' in the general case, it seems reasonable to conjecture that the the errors fall off as 
for large N. 

As a point of interest, it takes only 14 seconds to calculate 7 for a loop with A'' = A'^n + iVb = 32 on 
a Sun-4 workstation (SS2). The calculation time for 7 scales roughly as A'^^. The speed of this alogrithm 
makes it feasible to calculate 7 for loops with large numbers of segments A^a and Ni,. It is also possible to 
rapidly calculate 7 for very large numbers of loops with moderate values of Na and Nb- 

In this section we have shown that in all cases where previously published numerical methods have 
given reliable results for 7, these results are in good agreement with those given by our exact formulae. The 
large number of both piecewise linear and smooth loop trajectories for which our formulae have confirmed 
previously published results gives us confidence that our method is correct. In cases where our method yields 
result that disagree with previously published results, we have shown that our results are correct. We have 
shown this by identifying the errors in the previously published work and by showing that our results are 
consistent with those given by other independent methods such as the FFT method of Allen and Shellard 
[14] and/or a corrected implementation of the numerical method used by the original author(s). Since our 
results have been correct for every trajectory tested so far, we have confidence that our formulae provide a 
reliable method for calculating the power radiated in gravitational waves for arbitrary cosmic string loops. 

SECTION 8: CATALOG OF LOOPS 

This section gives a short catalog of piecewise linear loop trajectories and their 7 values. This catalog 
is intended to give a number of simple cases which might prove useful in testing future analytic or numerical 
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methods. The a- and b-loops which define these trajectories are regular polygons formed with small numbers 
of segments Na and A^^. With each pair of a- and b-loops we form a two parameter family of loop trajectories. 
The two parameters (</> and 9) describe the relative orientation of the a- and b-loop. There is nothing special 
about these loops other than that they are simple piecewise linear trajectories. 

The first set of trajectories we consider are defined by a- and b-loops consisting of 2 and 3 segments 
respectively. The three segment b-loop has the shape of an equilateral triangle. (The simplest case, where 
the a- and b-loops each have just two segments is discussed in the previous section.) The a-loop is taken 
to lie along the z-axis. One kink on the a-loop is positioned at the origin; the parameter u = at this 
kink. The other kink (at u = 1/2) is positioned above the first kink and has coordinates (0, 0, 1/2). For the 
b-loop, we again position one kink at the origin and set the parameter v = at that kink. The position of 
the other two kinks depends on the parameters ^ and 0. When (f) = = 0, the b-loop lies in the x-z plane. 
The kink a,t v = 1/3 has coordinates (—1/6, 0, y/3/6) and the kink at u = 2/3 has coordinates (1/6, 0, y/3/6). 
When ^ and 6 are not zero, the position of the b-loop is found as follows. First, the b-loop is rotated by the 
angle > about the z-axis (counter-clockwise when viewed from large positive z). After the (f) rotation, 
the b-loop is then rotated by the angle ^ > about the x-axis (counter-clockwise when viewed from large 
positive x). Values of 7 for the trajectories defined by these a- and b-loops are given in Table 1 for several 
values of the angles (p and 9. 



Table 1 



9\^ 


0° 


18° 


36° 


54° 


72° 


18° 


59.33 


59.80 


61.21 


63.51 


66.30 


36° 


53.81 


54.56 


56.86 


60.93 


67.45 


54° 


49.35 


50.15 


52.56 


56.40 


60.72 


72° 


46.70 


47.54 


50.12 


54.47 


60.36 


90° 


45.83 


46.70 


49.37 


54.12 


62.25 



In generating Table 1, we have been careful to avoid certain values of (j) and 9 for which the a- and b-loops 
have special relative positions. In particular, if the a- and b-loops are exactly co-planer, the operator D 
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becomes singular. In this case, accurate values of 7 may still be found by examining trajectories where the 
angles cj) and 9 deviate very slightly from their desired values. However, since the 7 values given in this 
section are meant to be "benchmark" values for future work, we have not included such cases here. In the 
following tables, certain pairs of angles are omitted for similar reasons. 

The second set of trajectories we consider are defined by a- and b-loops consisting of 3 segments each. 
Both the a- and b-loops are equilateral triangles. The position of the b-loop depends on the two parameters 
<j) and 9 in exactly the same way as the b-loop in the first set of trajectories. The a-loop is placed in the 
same position the b-loop has for parameter values cj) = 9 = Q. Values of 7 for the trajectories defined by 
these a- and b-loops are given in Table 2 for several values of the angles (j) and 9. 



Table 2 



e\4^ 


18° 


36° 


54° 


72° 


90° 


18° 


100.85 


90.65 


74.48 


64.80 


58.92 


36° 


82.00 


76.01 


70.97 


64.82 


59.98 


54° 


72.61 


66.51 


63.68 


61.54 


59.95 


72° 


67.04 


61.24 


59.12 


58.80 


60.02 


90° 


63.97 


58.53 


56.80 


57.12 


59.75 



The third set of trajectories we consider are defined by a- and b-loops consisting of 2 and 5 segments 
respectively. The two segment a-loop is identical to the a-loop used in the first set of trajectories. The 
b-loop is taken to be a pentagon. One kink on the pentagon is positioned at the origin and is chosen to 
have parameter value v = Q. When cj) = 9 = 0, the b-loop lies in the x-z plane, and is positioned so 
that the kink at t; = 1/5 has coordinates (— | cos(7r/5), 0, g sin(7r/5)), the kink at v = 2/5 has coordinates 
(i(sin(7r/10) — cos(7r/5)), 0, i(cos(7r/10) + sin(7r/5))), and so on. When (p and 6 are not equal to zero, the 
b-loop is rotated in exactly the same manner as for the previous two sets of loop trajectories. Values of 7 for 
the trajectories defined by these a- and b-loops are given in Table 3 for several values of the angles ^ and 9. 
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Table 3 



\ (p 


n° 
u 


1 

±c5 


oo 




79° 


18° 


63.52 


64.15 


66.04 


69.28 


74.52 


36° 


54.07 


54.99 


57.78 


62.31 


67.63 


54° 


47.69 


48.74 


52.02 


57.78 


66.36 


72° 


44.20 


45.30 


48.74 


54.98 


64.04 


90° 


43.10 


44.20 


47.69 


54.08 


63.95 



The fourth set of trajectories we consider are defined by a- and b-loops consisting of 5 and 3 segments 
respectively. The a-loop is a pentagon placed in the same position as the b-loop in the third trajectory set 
for parameter values 4> = 6 = Q. The b-loop is an equilateral triangle whose position is given in terms of 
the parameters (j) and 9 in exactly the same way as the b-loops used in the first and second trajectory sets. 
Values of 7 for the trajectories defined by these a- and b-loops are given in Table 4 for several values of the 
angles (j) and 9. 



Table 4 



9\cj> 


18° 


36° 


54° 


72° 


18° 


84.69 


75.43 


67.82 


62.44 


36° 


77.37 


71.03 


65.71 


61.86 


54° 


70.11 


65.69 


62.72 


60.42 


72° 


64.41 


61.37 


61.13 


61.33 


90° 


60.49 


58.21 


60.49 





The final set of trajectories we consider are defined by a- and b-loops consisting of 5 segments each. Both 
loops are taken to be pentagons. The a-loop is identical to the a-loop used in the fourth set of trajectories. 
The b-loop is in the same position as the a-loop when (f) = 9 = 0. When (f) and 9 are not zero, the b-loop is 
rotated in the same manner as in the previous sets of trajectories. Values of 7 for the trajectories defined 
by these a- and b-loops are given in Table 5 for several values of the angles (p and 9. 
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Table 5 



t> \ (p 




oo 




79° 


18° 


114.46 


94.04 


80.22 


68.84 


36° 


93.52 


82.49 


72.06 


65.40 


54° 


77.15 


72.94 


67.22 


62.74 


72° 


67.11 


64.76 


64.05 


62.98 


90° 


61.24 


59.47 


61.24 





The five sets of loop trajectories along with the 7 values given in this section are intended as "benchmark 
values" for future analytic or numerical work. 

CONCLUSION 

We have derived a new method for calculating the power emitted in gravitational radiation by cosmic 
string loops. This method yields an exact analytic formula in the case of piecewise linear cosmic string loops. 
By increasing the number of segments used, piecewise linear string loops can approximate any cosmic string 
loop arbitrarily closely. Our formula (derived in sections 5 and 6) involves nothing more complicated than 
log and arctangent functions. No numerical integrations are required. Further, since our formula is exact, 
there is no need to estimate any contribution to 7 from the "tail" of an infinite sum. The error introduced 
when approximating smooth loop trajectories by piecewise linear trajectories has been investigated. It is 
found that this error typically falls off as for large A^, although in at least some cases it falls off faster, 
as N~^. We believe that for "generic" loops the error scales as N~^. Using a computer to evaluate the 
approximately N'^ terms in our formula, we can determine values of 7 more accurately and more efficiently 
than by previously published methods. 

We have tested the results of our formula against all previously published radiation rates for different 
loop trajectories. Section 7 contains a detailed comparison of the results given by our new method to those 
reported by previous authors. In every case, our formula is found to give the correct result. In many cases our 
results are in good agreement with the published data. However there are also a number of cases where our 
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results do not agree with those previously published. There are, in fact, a number of cases where conflicting 
results have been published for the same trajectories. In the cases where a disagreement was found, we have 
identified the errors made in the published work which led to the incorrect results. In most cases, the errors 
in the published values of 7 are a result of underestimating the contribution of the tail of the infinite sum in 
(7.3). The incorrect values of 7 which have been published are typically 25 to 50 percent below the correct 
results. We are confident that our formula gives the correct results because, in every case, we have shown 
our results to be consistent with those given by independent methods. 

We intend to use this exact formula in future work, for example to repeat some of the work of Scherrer, 
Quashnock, Spergel and Press [12] concerning the distribution of 7 values of non-self-intersecting loops. In 
addition, we plan to show how this formula may be modified to yield similar analytic results for the linear 
momentum radiated by cosmic string loops [19]. 
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FIGURE CAPTIONS 



Figure 1. For piecewise linear loops, the a- and b- loops consist of straight segments. The segments are 
joined together at kinks where a'(u) and b'{v) are discontinuous. The kinks on the a-loop are labeled by 
the index i. The spatial position of kink i is a,. The value of the coordinate u at the ith kink is The 
segments on the a-loop are also labeled by the index i; the ith segment being the one between the ith and 
the {i + l)th kink. 

Figure 2. The delta function, 5{f{x,y,z)), which appears in (5.2) only has support when f{x,y,z) = 0. 
Solving f{x, y,z) =0 for z{x, y) we find that the surface z{x, y) consists of a pair of disconnected hyperbolic 
sheets. The hyperbolic sheets are separated by the plane (iix + /32?y + /34 = 0, where the denominator of (6.4) 
vanishes. The intersection of a ^ =constant plane with these sheets will be a hyperbola in the x-y plane. 
The integration volume for (5.2) is a box with opposite corners (0, 0, 0) and (Aa;, Ay, A2;). 

Figure 3. The integrals (6.9) in the expression for the radiated power are over the rectangular region bounded 
by the corners (0,0) and (Aa;, At/). Each integral contains three step functions. In the region of integration, 
the first two step functions are both non-zero only in the shaded regions between the "Top" and "Bottom" 
hyperbolae. The third step function is on in the region which includes the "future" branches if s = 1 or the 
"past" branches if s = — 1. The vertical asymptote of the Bottom (Top) hyperbola is shown as a dashed line 
that lies at x = (ps {x = 

Figure 4. The five diflPerent fundamental types of integrals that arise when evaluating (6.9). The type depends 
upon the relative positions of the two hyperbolic curves y{x,0) (labeled B for "Bottom") and y{x,Az) 
(labeled T for "Top") and the rectangular region bounded by the opposite corners (0, 0) and (Ax, Ay). The 
five possibilities are shown using the future branches of the hyperbolae (the case for s = -|-1). Each type of 
integral may also occur with the past branches of the hyperbolae (for s = —1). 

Figure 5. The solid line is a plot of the analytic formula (7.1) for 7((/>) for piecewise linear trajectories in 
which the a- and b-loops are composed of just two segments each. The dots show values of 7 given by our 
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formulae for a range of angles </>. The agreement between the two results is so close that the dots had to be 
enlarged in the figure to distinguish them from the solid line plot of (7.1). 

Figure 6. The solid curves show numerical values of 7 for piecewise linear loop approximations of M = 1, 
L = 3,5 Burden loops. The piecewise linear a-loops were constructed to have Na = 16L segments. The 
piecewise linear b- loops each had Ni, = 16M =16 segments. The open circles show the published results of 
Burden [10], the crosses show the published results of Quashnock and Spergel [13], and the triangles show 
the results of the FFT method of Allen and Shellard [14]. We find excellent agreement among all four sets 
of results. 

Figure 7. Numerical values of 7 for some L = M = 1 Burden loops. The solid line shows the results of our 
new method. The open circles show the results of Burden, the crosses show the results of Vachaspati and 
Vilenkin, and the triangles show the results of the FFT method. There is excellent agreement among these 
four sets of results. Durrer's results are shown as open diamonds. 

Figure 8. The 7 values reported by Durrer (open diamonds) compared to the sum of the first 50 terms (small 
crosses) in (7.3) for several L = M = 1 Burden loops. Including an estimate of the contribution to 7 from 
the infinite tail of the sum results in significantly larger values of 7 (large crosses). The solid line shows the 
values of 7 found by our formulae. 

Figure 9. Values of 7 for piecewise linear approximations to the L = M = 1 Burden loops using different 
numbers of segments {Na, Nh). The number of segments {Na, A^fe) used for the curves A, B, C, D are 
respectively, (6,5), (11,10), (16,15) and (36,35). The Burden loops are accurately approximated over the 
range (j) e [10, 160] degrees when Na > 16 and N^ > 15. Regions ^ < 10 and ^ > 160 where 7 begins to 
diverge require a larger number of segments before the approximation becomes accurate. 

Figure 10. Relative errors in 7 for the {Na, Ni,)={6,5), (11,10), (16,15) and (26,25) piecewise linear loop 
approximations of the L = M = 1 Burden loops with respect to the {Na, Ni,)={36,35) approximation. 
The relative error e of each set of loops is calculated by summing [(7^ — 7|')/(7^ + 7|')| over the values 
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(j) = 5,10, ... , 175 degrees and then dividing by the number of terms in the sum. Here, a denotes which set 
of loops are being compared (i.e. (A^o, iVb)=(6,5), (11,10), etc.) and E denotes the {Na, A^b)=(36,35) loops. 
Most of the error is due to loops with values of (p near or 180 degrees. Increasing the number of segments 
from {Na, Nh)={16,15) to (36,35) causes the average value of 7 to change by less than 3 percent. 

Figure 11. Values of 7 are shown for the loop trajectories given in (7.4) with a = 0.5. The results given by 
our formulae are shown by the sohd hne. The results of the FFT method are shown as triangles. There is 
good agreement between the FFT method and our new method. Durrer's results are shown as dots while 
the results of Vachaspati and Vilenkin are shown as crosses. These last two sets of results are inaccurate 
because the rate of convergence of the sum (7.3) was estimated incorrectly. The open circles show the sum 
of the first 300 terms of (7.3) as given by the FFT method. These circles should be taken as lower bounds 
on the values of 7. 

Figure 12. The power spectrum for the trajectory (7.4) with a = 0.5 and 4> ~ 11/2 found using the FFT 
method. The P„ are shown in units of GjjL^. The slope of the curve for large n shows that the sum (7.3) falls 
off as the power law P„ a n~^'^^ for large n. The logarithm is to base 10. 

Figure 13. Values of 7 are compared for the loop trajectories given in (7.4) with a = 0.01. The results of our 
formulae are shown by the solid line. The results of the traditional numerical method are shown as crosses. 
The two methods are in good agreement. The values of 7 when a = are shown by the dashed line. 

Figure 14. Values of 7 given by our formulae for the trajectories (7.5) for three sets of parameters (p, q) 
and a range of angles 4>. Curves A, B and C give the results for trajectories with {p,q) =(0.6,0.4), (0.4,0.8) 
and (0.9,0.9) respectively. These results are consistent with the claim by Garfinkle and Vachaspati that the 
trajectories (7.5) have 7 values on the order of 100. The results of the FFT method with (p, g) =(0.9,0.9) 
are shown by the triangles, and should be compared to curve "C". There is good agreement between the 
two methods. 

Figure 15. The 7 values from a series of increasingly accurate piecewise linear loop approximations to two 
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[L = M = 1) Burden loops. The total number of linear segments in each approximation is given by N. Both 
sets of 7 values quickly converge to their asymptotic limits. 

Figure 16. The function A(Af) for two (L = M = 1) Burden loops. Each point shows the difference between 
the 7 value for a piecewise linear approximation with A'' segments and the numerical estimate of 7 in the 
N = 00 limit. Both of the solid lines have a slope of -1, showing that the errors fall off as N~^. The 
logarithm is to base 10. 
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